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ABSTRACT 



Using a large set of XMM-Newton observations we searched for long term spectral and flux variability of the isolated neutron star RX 
JI856.5-3754 in the time interval from April 2002 to October 2011. This is the brightest and most extensively observed source of a 
small group of nearby, thermally emitting isolated neutron stars, of which at least one member (RX J0720.4-3 125, see e.g. Hohle et al., 
2010) has shown long term variability. A detailed analysis of the data obtained with the EPIC-pn camera in the 0.15 - 1.2 keV energy 
range reveals only small variations in the temperature derived with a single blackbody fit (of the order of ~ 1% around an average 
value of fcr°° ~ 61 eV). Such variations appear to be correlated with the position of the source on the detector and can be ascribed to 
an instrumental effect, most likely a spatial dependence of the channel to energy relation. For the sampled instrumental coordinates, 
we quantify this effect as variations of ~ 4% and ~ 15 eV in the gain slope and offset, respectively. Selecting only a homogeneous 
subset of observations, with the source imaged at the same detector position, we find no evidence for spectral or flux variations of 
RX JI856.5-3754 from Mai'ch 2005 to present-day with limits of AfcT" < 0.5% and A/x < 3% (0.15 - 1.2 keV range), with 3cr 
confidence. A slightly higher temperature (kT'" ~ 61.5 eV, compared to the average value kT°° ~ 61 eV) was instead measured in 

April 2002. If this difference is not of instrumental origin, it implies a rate of variation 0. 15 eV yr"' between April 2002 and March 

2005. The high-statistics source spectrum from the selected observations is best fitted with the sum of two blackbody models, with 
temperatures kT^ = 62.4+jJ-6 eV and kT^ = 38.9+^-^ 

narrow or broad spectral features are detected, with upper limits of ~ 6 eV on their equivalent width. 
Key words. Stars: neutron. Stars: Individual: RX J1856.5-3754 



38.9^to eV, which can also account for the flux seen in the optical band. No significant 



1. INTRODUCTION 

RX J1856.5-3754 (J1856 hereafter) is the prototype of the so- 
called XDINSs, i.e. the seven X-ray Dim Isolated Neutron S tars 
discovered by the ROSAT satellite (e.g. IWalter et al-lligg e*. see 
also .Haberl 2007; Tur olla 2009 for reviews). These are nearby, 
d < 300 pc, radio-quiet isolated neutron stars characterized by 
thermal spectra, with temperature^ ~ 50 - 100 eV and lu- 
minosities Lx ~ lO-'' - lO-'^ erg s"' . Six out of the seven XDINSs 
exhibit broad spectral features at energies between 270 and 700 
eV, with equivalent widths of several tens of eV, interpreted as 
proton cyclotron lines or atomic transitions in strong magnetic 
fields, B ~ 10'^ G. Spin periods of the XDINSs are in the 3-12 
s range and, assuming magneto-dipole braking, their measured 
spin-down rates, of ~ lO"''* - 10"'^ s s"', imply magnetic fields 
B ~ 10'^ - 10'"* G, in broad agreement with those inferred from 
the spectral features. All XDINSs have very faint optical/UV 
counterparts, my ~ 26 - 27: the optical/UV flux exceeds that 
expected fro m the Rayleigh-Jea ns tail of the X-ray blackbody 
(BB, see e.g. lKaplan etairi201 li) . 

Intriguingly, the second brightest XDINS, RX J0720.4-3125 
(J0720 hereafter) showed significant long term variations of its 



measured by an observer at infinity. 



spectral properties (iHohle et al.L 1201 Ol and references therein). 
The nature of these changes is still under deb ate and can possibly 
be related to a precession of its spin axis (e.g.'Haberl et al.V200^ 
or to a glitch-like episode (van Kerkwiik et al., 2007). Another 
interesting fact is that XDINSs share some of their properties 
with magnetars. In particular, their spin periods fall in the same 
range, while the magnetic fields of XDINSs are intermediate be- 
tween those of magnetars, B ~ 10'"*- 10'^ G, and normal pulsars, 
B ~ 10" - 10'^ G. Also, XDINSs seem to be older than magne- 
tars, with spindown ages of ~ 10^-10^ years against ~ 10^-10^ 
years (e.g.'M ereghettiLl2008h . These similarities suggest a possi- 
ble evolutionary Unk between the two classes of isolated neutron 
stars (INSs). 

J 1856 has the lowe st pulsed fraction, ~ 1% 
(iTiengo & Mereghetti'. '2007), and the weakest magnetic field, 
B ~ 1.5 X 10"^ G (van Kerkwijk & Kaplan, 2008) among known 
XDINSs. Being the brightest, fx ~ 1.5 x 10 " erg_s_^_crn_^ 
and closest of the category, d = 123![^pc (e.g. ' Walter et"an. 
2010), it stands among the best candidates to probe the mass 
and radius of a neutron star, which in turn would be of 
paramount importance to constrain the equation of state (EOS) 
of matter at nuclear densities. Analyzing the data from a 57 
ks XMM-Newton observation and a ~ 500 ks observation 
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by Chandra, iBurwitz et alJ (l2003h suggested that atmosphere 
models with heavy elements can be ruled out, given the lack 
of any spectral features. Also, they ruled out non-magnetic 
and fully ionized hydrogen atmosphere models, since both 
would over-predict the optical flux. On the other hand, they 
found that the broadband (optical + X-ray) spectrum is well 
fitted by a two component BB, where the hotter component, 
kT^ ^ 63.5eV(/?^ ^ 4.4(t//120pc)km), is responsible for 
the X-ray emission while the cold component, kTopi < 33 eV 
(R'^p, > 17 (t//120pc)km), accounts only for the optical emis- 
sion. Alt erna t ive models, like emi s sion f rom condensed matter 
surface dLail 120011; |Turolla et al.', '2004), or partially ionized 
hydrogen atmospheres (Ho et al ., 2007) cannot be ruled out. 

Thanks to its bright and presumably steady emission, J 1856 
has been targeted routinely for calibration purposes by XMM- 
Newton in the last decade. The large amount of collected data 
allows us to precisely characterize its spectral evolution on time 
scales from months to ~ 10 years. At the same time, the stabil- 
ity of the detectors on board XMM-Newton in the same time- 
frame can be scrutinized as well, as done with the data taken 
before 2006 in Haberl (2007). Here we report on the analysis 
of all XMM-Newton observations of J 1856 pe rformed so fa r in 
imaging mode with the EPIC-pn camera ( St riider et a 1I I2OOI on 
board XMM-Newton. The reduction and the analysis of the data 
is described in Section |2] Results are also illustrated in Section 
|2l We constrain the detector stability and obtain an upper limit 
for the spectral variations of the star and compare it with those 
of other isolated neutron stars reported in the literature. A dis- 
cussion on the astrophysical implications of these results is also 
given (Section[3]|. 

2. DATA REDUCTION AND SPECTRAL ANALYSIS 

We use all imaging data collected with the EPIC-pn camera from 
April 2002 to October 20 11. Data from the two EPIC-MOS de- 
tectors are not considered in this analysis because the MOS ef- 
fective area at soft X-ray energies is much smaller than that of 
the pn and the MOS cameras are known to be less stable on 
long term (Read et al., 2006). Table [1] reports the date, identifi- 
cation number, net exposure time, and number of background- 
subtracted counts of J 1856 for each observation. All observa- 
tions were performed in Small Window mode with the thin filter. 
Raw data are processed with the XMM-Newton science analysis 
package (SAS) vl 1.0, using the epproc pipeline. We then build 
light curves above 10 keV, with bin size of 100 s, in order to 
identify soft proton flares. Based on these light curves, we clean 
the data by removing time intervals with count rate higher than 
4cr above the mean rate. Source spectra are then extracted from 
a circular region of 30" radius, selecting only single pixel events 
(i.e. PATTERN=0 in the xmmselect task), in order to obtain 
event lists of the highest possible quality. The extracted spec- 
tra are binned in order to have at least 30 cou nts per bin. Th e 
spectral analysis is performed with XSPEC 12 (lArnaudLll996h . 
selecting photon energies in the 0.15 - 1.2 keV range. 

We perform a simultaneous fit of all the 20 spectra with 
a single absorbed BB, using a phabs'-bbodyrad model and 
adopting the abundances of Anders & Grevesse ( 1989). The true 
spectra of XDINSs can be more complex than a single BB. 
However, we are interested in relative variations in the emis- 
sion of J 1856. Thus, we rely on the parameters of the single 
BB in order to quantify these variations. The best fit gives a 
BB temperature kT'" = 61.30 + 0.04 eV, and a column density 
A^H = (5.84 + 0.04) X 10''^ cm-2 (j^l = 1.78 with 2784 degrees of 
freedom). The residuals show large systematic deviations, espe- 



cially above ~ 0.5 keV (Fig.lT] upper panel), which change from 
observation to observation, indicating that some of the spectral 
parameters vary between different observations. Therefore, we 
repeat the spectral fit leaving the BB temperature and normal- 
ization as free parameters for each observation. In this case we 
obtain a better xl = 1-37 with 2746 degrees of freedom. The 
resulting parameters are shown in Fig.|2](panels a, b, c). The fit 
of the 20 spectra with the column density as free parameter, and 
a unique temperature and normalization for all observations, is 
not as good as the previous one, xl - 1-58 with 2765 degrees 
of freedom. The resulting column densities are shown in Fig. |2] 
(panel d). The count rate shows variations of ~ 1-2 percent 
around a mean of ~ 7.46 cts s ', but with no evidence of a de- 
fined trend and without correlation with the BB temperature or 
normalization. On the other hand, the hydrogen column density 
is clearly anti-correlated with the count rate (compare panels a 
and d in Fig.|2]l, indicating that the changes in the count rate are 
mainly determined by diff'erences in the so ftest spectral band. As 
akeady noticed bv lStuhlinger et all (1201 Ol) . the temperature also 
shows variations of ~ 1-2 percent (~ 1.5 - 2eV), which are 
anti-correlated with the normalization. We note that these varia- 
tions have a seasonal pattern with a period of ~ 12 months, sug- 
gesting that they are likely related to instrumental eff'ects rather 
than being intrinsic changes of the source emission. To verify 
this hypothesis we study the relation of the BB temperature with 
the position of the source centroid on the detector (RAWX and 
RAWY coordinates). With the exception of the three off'-axis ob- 
servations (observations H, P and Q, which includes the four 
separate pointings performed in October 2011) and one in an 
intermediate position (observation F), in most of the observa- 
tions the source centroid lies in two close but distinct regions 
of the detector (Fig.[3]l, one at RAWX~36-37, RAWY~192, and 
the other one at RAWX~38, RAWY~190 (hereafter 'soft' and 
'hard' region, respectively, according to the relative hardness of 
their spectra, as outlined below). This is due to the diff'erent roll 
angle of the satellite in the Spring and Fall visibility intervals for 
the J 1856 sky region. 

The best fit BB temperature is systematically higher for ob- 
servations with the source in the 'hard' region (upper right and 
lower left panels of Fig. O, and in the 3 most extreme ofi" axis ob- 
servations at RAWY=168, RAWX=16 and RAWX=56. In the pn 
camera the signal is read independently from each column along 
the RAWY position. This implies that small inaccuracies in ei- 
ther the gain or the charge transfer inefficiency (CTI) correction 
along the readout direction could be the underlying cause for the 
observed spectral differences. Thus we fit again all the spectra 
with a unique BB model, but this time allowing a variation of 
the gain parameters. These parameters, slope and offset, define 
the relation between the energy scale in the response matrix of 
each observation with respect to that of an observation taken as 
reference: E'- - E/slopei - offsetj. We assume the spectrum of 
the first observation (observation A) as reference, i.e. we fix its 
slope and offset at 1 and 0, respectively. In this case we obtain 
xl - ^ -32 with 2746 degrees of freedom, with a net reduction of 
the scatter of fit residuals for all spectra (Fig. [1] lower panel) and 
values of the gain parameters ranging from ~ 0.986 to ~ 1.028 
for the slope and from ~ -6.6 eV to ~ 8.1 eV for the offset. This 
implies that there could be variations up to ~ 4% of the slope and 
of ~ 15 eV of the offset throughout the -10 years and source 
positions sampled by the XMM-Newton observations (see Fig. 



' Note that the PN channel energy has a width of ~ 5 eV, which 
means that the photon energy as detected by the instrument is not known 
better than that. 
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Table 1. Log of all EPIC-pn observations of RX J1856.5-3754 performed in Small Window Mode. 



Observation 


Date 


Obs.ID 


Net Exposure Time 
[s] 


Counts 


A 


2002-04-08 


0106260101 


40030 


302601 


B 


2004-09-24 


0165971601 


22960 


174955 


C 


2005-03-23 


0165971901 


14190 


106949 


D 


2005-09-24 


0165972001 


22700 


168609 


E 


2006-03-26 


0165972101 


48340 


361845 


F 


2006-10-24 


0412600101 


49820 


374018 


G 


2007-03-14 


0412600201 


26940 


204513 


H 


2007-03-25 


0415180101 


16540 


126883 


I 


2007-10-04 


0412600301 


23920 


179449 


J 


2008-03-13 


0412600401 


33080 


245923 


K 


2008-10-04 


0412600601 


43580 


330553 


L 


2009-03-19 


0412600701 


47510 


353698 


M 


2009-10-07 


0412600801 


40780 


306449 


N 


2010-03-22 


0412600901 


48380 


360306 





2010-09-29 


0412601101 


47870 


356580 


P 


2011-03-14 


0412601301 


45210 


342204 


Qi 


2011-10-05* 


0412601501 


17490 


134747 


Q2 






16310 


124431 


Q2 






16800 


127176 


Q4 






17920 


128549 



(*) The observation of October 201 1 has been divided in 4 pointings, 



with the source at different positions on the detector. 



[3]l. This is a conservative estimate since we considered the two 
most extreme spectra and assumed no intrinsic spectral variabil- 
ity in J 1856. 

2.1. CONSTRAINTS ON INTRINSIC SPECTRAL 
VARIATIONS 

To investigate possible variations of the spectral properties of 
J 1856 during the nine year-long time-span of the monitoring 
campaign, we have to compare only homogeneous data sets, i.e. 
when the source was located at approximately the same detector 
coordinates, in order to reduce the incidence of the systematic 
uncertainty reported in the previous section. 

First, we study temperature variations considering the 8 
spectra of the observations with the source in the 'soft' region. 
Fitting a linear function to the derived temperatures, we obtain a 

negative (cooling) rate of variation of 0.035 + 0.010 eV yr"' 

(Fig. |4]l. However, the fit is unacceptable (xl = 5.16 for 6 de- 
grees of freedom) because of the first data point (observation A), 
which differs significantly from the others. Assuming no varia- 
tions of the instrument response between April 2002 and March 
2005, the 'anomalous' temperature inferred from the older spec- 
trum may be explained by an intrinsic change of the source emis- 
sion, like those observed in J0720, implying a relatively rapid 
drop in the temperature, of ~ 0.5 eV in less than 3 years. In alter- 
native, some subtle variations of the instrument (energy response 
of the pixel column) may have occurred in the same time-frame. 
Excluding the first data point we obtain a better fit (xl - 1 .74 
for 5 degrees of freedom) and the rate of temperature variations 
of ~ 0.023 ± 0.015 eVyr"', which is de facto consistent with a 
constant temperature at 2cr level (Fig. [Hi. 

If instead we use the data from the 5 'hard' observa- 
tions, we obtain a rate of temperature variation of ~ 0.044 + 
0.026 eVyr"', again consistent with a constant temperature at 
2cr level. However the value obtained is affected by larger un- 
certainty since, as it can be seen from the middle panel of Fig. 
121 in this case there is a larger scatter of the kT°° values with re- 



spect to those found for the 'soft' observations. The origin of this 
larger scatter is unclear and could be possibly related to a gain 
instability of the pixel columns within the source point spread 
function. In any case, owing to their homogeneity and longer to- 
tal exposure, throughout the rest of the paper we will make use 
of only the 'soft' data in order to constrain the spectral properties 
of J1856. 

2.2. BEYOND THE SINGLE BLACKBODY 

In order to obtain a single spectrum with high statistics, we 
merge together the 7 spectra of the 'soft' group (Fig. [3]) ex- 
cluding, as before, the anomalous observation of April 2002 
(observation A). The resulting spectrum contains ~ 1.9 x 10^ 
background-subtracted counts and coiTesponds to a net integra- 
tion time of ~ 254 kiloseconds. 

We first try fits with a single BB (Fig.|5] upper panel). In or- 
der to obtain a formally acceptable fit (null hypothesis probabil- 
ity > 0.1) we add a systematic eiTor of 1 .5% (systematic param- 
eter in XSPEC). This is an energy-independent systematic error 
added to the model in XSPEC and accounts for the uncertainties 
on the spectral model and for likely residual calibration inaccu- 
racies. We obtain a reduced = 1.12 for 176 degrees of free- 
dom. The resulting temperature is kT°° = 61.5 + 0.1 eV, to which 
coiTesponds an emission radius = 5.0 + 0.1 km for a distance 
of 120 pc. The column density is A^h = (4.8 + 0.2) x lO'*^ cm^^. 
The fit residuals exhibit systematic deviations both at the lowest 
and highest ends of the energy range (Fig.|5] upper panel), which 
may be indicative of the inability of the single temperature BB 
to adequately fit the spectrum of J 1856. 

The presence of a second BB at diff'erent temperature is ex - 
pected to explain pulsations (e.g. iTiengo & MereghettiL l2007h . 
We therefore fit the merged spectrum with a two BB model 
(Fig. |5] lower panel). In this case we obtain an acceptable fit 
by adding a systematic error of 0.6%, which returns a;i'^ = 1.11. 
The column density is now A^h - (12.9 + 2.2) x 10'^ cm"^. The 
hard BB has a temperature of kT^ = 62.4;^Q4eV with emis- 
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sion radius of R'^ = 4.7+°-2 (c//120pc) km, while the soft BB 
has a temperature kT^ - 38.9;^2 9 ™d emission radius Rf - 
11.8;^Q4((i/120pc)km. We summarize these resuhs in Table|2l 
Interestingly, the contribution of both the soft and hard X-ray 
blackbodies to the optical flux can account for the exc ess ob- 
serve d in J 1856 (van Kerkwiik & Kulkarni, 2001; Kapl an et al.L 
I2OI ih . Considering the best-fit parameters for the double BB 
model, the contribution of the soft BB at optical wavelengths 
is ~ 4 times larger than that of the hard BB. This implies an 
overall increase of a factor ~ 5 with respect to the optical flux 
expected from the Rayleigh- Jeans tail of the hard BB alone (Fig. 
|6l). This value is consistent, within uncertainties, to the factor 
~ 7 obtained from the comparison between opt ical/UV photom- 
etry and the extrapo lation of the X-ray data ( Bur witz et aUl2003t 
lKaplanetalll2011h . 

We note three dips at ~ 0.3, 0.4 and 0.6 keV in the residuals 
of both single and two-component BB fits. Since their width is 
smaller than the energy resolution of the pn detector in the range 
of interest, they are most likely due t o non perfect i nstrumental 
calibrations, as already suggested bv iHaberll (l2007h . We search 
for other possible features in the spectrum of J 1 856, adopting the 
single BB model with no systematic errors in order to facilitate 
the comparison with the spectral features found in other XDINS. 
We investigate the presence of narrow absorption features from 
300 to 700 eV at 100 eV intervals. We use the gaussian model 
from the XSPEC library, imposing fixed line energy and widtlfl 
(cr = in XSPEC) and negative normalization. We find a feature 
at 400 eV, most likely one of the aforementioned dips and there- 
fore of instrumental origin, with equivalent width of 2.0 + 0.5 eV 
at 4cr confidence. No other narrow features are found at the same 
confidence level, corresponding to a maximum equivalent width 
of ~ 6eV. Typically, XDINSs exhibit broad absorption lines. 
Hence, we look for broad features (fixed line width cr = 0. 1 ke V) 
but with line energy free to vary. We find a feature at 370 + 15 eV 
and equivalent width 9 + 2 eV at Act confidence. However, the re- 
ality of this line is debatable and could arise from an attempt to 
fit a non purely Planckian spectral continuum. In fact, if we adopt 
the two-component BB as starting model, no broad features are 
found at 4cr confidence level. 



3. DISCUSSION 

A spectral analysis of all the imaging observations of the isolated 
neutron star RX J1856.5-3754, performed from 2002 to 2011 
with the EPIC-pn instrument on-board the XMM-Newton satel- 
lite, revealed small amplitude variations, of the order of ~ 1 -2%, 
quantified by changes in the fit parameters obtained with a sin- 
gle absorbed BB model. The long term trend of these variations 
suggested an instrumental origin, since there is a correlation of 
the fit parameters with the position of the source image on the 
detector. This is likely related to a non-uniform energy response 
between different channels of the readout electronics. The re- 
lated uncertainty has been quantified with a gain fit in XSPEC, 
resulting in variations of the gain slope and offset of ~ 4% and 
~ 15 eV, respectively, over the pn detector positions covered by 
the J 1856 observations (Fig. [3]). 

Using only homogeneous data, we found that the rate of tem- 
perature variations, in the framework of the single BB model, is 
compatible with a constant at 2cr confidence level over the last 



~5 years of the monitoring campaign. A higher temperature was 
recorded on April 2002. If not caused by subtle alterations of 
the instrument response, this difference would imply that also 
J 1856 undergoes spectral changes, albeit the magnitude of these 
changes is much smaller than that observed in J0720. The ob- 
served temperature changed of ~ 0.5 eV in three years (or less, 
considering the lack of coverage between April 2002 and March 

2005) corresponds to a rate of 0.15 eV yr"'. In any case, the 

long term behavior of J 1856 is markedly different from that of 
J0720. The latter underwent substantial and continuous changes 
in its observed spectral properties during the years. Apart the 
slightly higher temperature recorded in April 2002, J 1856 ex- 
hibits instead a steady behavior, with no spectral variations in 
the last 5 years. 

By merging all the homogeneous data-sets we obtained a 
high statistics spectrum with which we tested the validity of sin- 
gle and two BB models and checked for the presence of features 
in the spectrum of J 1856. Apart some narrow features likely 
due to calibra tion issues and already reported in the literature 
(IHaberll l2007l) . we found no convincing evidence for broad or 
narrow absorption features. The absence of lines thus remains 
one of the distinctive p roperties of J 1856 among XDINSs (see 
also lHohle et al.Ll201 ih . The two BB model returns a better fit to 
the data and is more physically justified by the observed pulsa- 
tions. The presence of a second, cooler BB component was in- 
voked by Pons et al. (2002) to explain the optical/UV data but its 
contribution to the X-ray flux was considered negligible because 
of the apparent lack of pul sations reported for J 1 856 at that time 
(Braie & Romani, 200^ iBurwitz et all l2003t iTriimper et all 
2004). This implied an upper limit on the temperature of the 
cold BB, kT^p, < 33 eV, which in turn yielded a lower limit on 

the stellar radius, R'^ = ((Rop,f + ( R^)^y'^ > 16(J/120pc) km. 
An even larger value was found by Hamb arvan et al.l (1201 1) for 
RBS 1223 (IRXS J130848.6-I-212708), another XDINS. It must 
be noted however that they assumed a condensed surface with a 
thin hydrogen atmosphere model to fit the spectrum of this neu- 
tron star. Thus, the measure of the radius resulting from the fit is 
larger than that obtained with a BB spectrum. From our analysis 
we found instead a smaller radius, 7?°° = {{RTf + {Rl^f y'^ = 
12.7'*'n2 (f//120p c)km, consiste nt at Icr level with the lower 
limit reported bv ITriimper et al.l Thus, our result allows a wider 



Table 2. Summary of fit parameters for the joint spectrum. 



Parameter Single BB Two BB 



^ The condition cr = implies that the line width is smaller than the 
energy resolution. The width a is thus automatically adjusted to match 
the minimum width imposed by the instrument resolution. 



Nh 
[lO'^cm-^] 

[eV] 

[km] 

kTT 
[eV] 

R7 
[km] 



4.8!«;2 12.9!|2 



615+"' 62 4+"'' 



c rv+0.1 47+0.2 
-^•"-0.1 -0.3 



38.9!^;^ 



1.5% 0.6% 
1.12 1.11 
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range of radii and is not much constraining for the EOS. In any 
case, we stress that the inferred parameters of a second BB in the 
X-ray spectrum of J 1856 must be taken with caution. In addition 
to the statistical and systematic errors reported above, the nor- 
malization and thus the radius of soft BB component is affected 
by systematic uncertainties because of the not well constrained 
energy redistribution below ~ 0.4 keV. Large systematic devia- 
tions are expected also because J 1856 is one of the main calibra- 
tion targets used to determine the instrumental energy response 
and so the pn response matrix is currently tuned to reproduce 
a X-ray spectrum of J 1856 not in contrast with its optical/UV 
flux, if a double BB models is assumed. With the caveat of the 
aforementioned uncertainties, the extrapolation of the cold BB 
acco unts also for a large fraction of the reported optical excess 
(e.g. iKaplan et al.L |20I I h. avoiding the need for more complex 
models proposed to explain the observed optical/UV flux. 



4. CONCLUSIONS 

The results presented in this work show that the small amplitude 
variations in the spectral parameters of the isolated neutron star 
RX J 1856.5-3754, obtained by fitting its spectrum with a single 
absorbed BB model, are due to a non-uniform energy response of 
the EPIC-pn camera. Once this instrumental effect is taken into 
account, the upper limits on the relative temperature and flux 
variations in the period March 2005-present are AA:r°° < 0.5% 
and Afx < 3%, respectively. These can be taken as a measure of 
the source H-instrument long term stability. A higher temperature 
was instead observed in April 2002. If due to an intrinsic change 
of the source spectral properties, it would imply that variations 
on time-scales of years, like those observed in the other XDINS 
RX J0720.4-3125, might be a common feature of this class of 
sources. 

We point out that, being the two brightest and most observed 
XDINSs, a continuous monitoring of these two sources will help 
to characterize the differences and analogies in their long term 
behavior and thus will give insights about the physical condi- 
tions on their surfaces, like e.g. the magnetic field and tempera- 
ture distributions. In the case of J 1856, the steady emission al- 
lowed us to sum all the data taken with the source in the same 
position on the detector, in order to obtain a spectrum with high 
count statistics. The resulting spectrum is best fitted by a two 
BB model, which is also more justified than a single BB by the 
observation of pulsations at X-ray energies. Also, the extrapo- 
lation at optical wavelengths of the emission from the two BBs 
accounts for the excess reported in the literature. However, due 
to uncertainties in the calibration of the pn at low energies, the 
radius of the star is not well constrained. 

Finally, the results presented here suggest that J 1856 is possi- 
bly the best choice as calibration source for instruments observ- 
ing in the soft X-ray domain. Moreover, the EPIC-pn detector on 
board XMM-Newton is extremely stable over a long time-span 
(see above) and so, after taking into account a subtle spatial de- 
pendence of its spectral response, it can be considered as a ref- 
erence instrument for the study of long term spectral variability 
of X-ray sources. 
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Fig. 2. Long term evolution of the spectral parameters obtained with a single BB fit in the 0.15 - 1.2 keV range. Panels a, b, 
and c, show the count rate, temperature, and BB normalization of a simultaneous fit with the column density fixed to a common 
value. Panel d indicates the column density obtained when, instead, the BB parameters are fixed to a common value for all the 
observations (see text). Diamonds, triangles and squares represent 'soft', 'hard' and ofF-axis observations, respectively. The plus 
symbol is the October 2006 observation, which is at an intermediate detector position (see Fig. O. Error bars correspond to 3cr 
confidence intervals. 
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Fig. 3. Lower right panel: position of the source in detector coordinates (RAWX,RAWY). The region in the dashed square is 
enlarged in the upper left panel. The symbols indicate the different observations: 'soft' region (diamonds), 'hard' region (triangles), 
observation F (plus), off-axis observations (squares). The other two panels indicate the temperature as a function of the source 
position. 
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Fig. 5. (Top) Single BB fit of the merged spectrum obtained from observations of the 'soft' group. (Bottom) Same as Top panel but 
with a two BB model. 
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Fig. 6. Broadband spectrum of J1856. The dashed line is the extrapolation of the (unabsorbed) X-ray BB from'Burwitz et al.' ("20031). 
The dotted lines indicate the two BB components obtained from the sum of all the homogeneous observations with the source in the 
'soft' region, and the solid line is their sum. The shaded area marks the Icr confid ence region for the best fit model. The triangles 
represent the optical/UV data obtained from HST photometry (iKaplan et al.L 120111) . 
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